Observation of Density Segregation inside Migrating Dunes 
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Spatiotemporal patterns in nature, such as ripples or dunes, formed by a fluid streaming over a 
sandy surface show complex behavior despite their simple forms. Below the surface, the granular 
structure of the sand particles is subject to self-organization processes, exhibiting such phenomena 
as reverse grading when larger particles are found on top of smaller ones. Here we report results 
of an experimental investigation with downscaled model dunes revealing that, if the particles differ 
not in size but in density, the heavier particles, surprisingly, accumulate in the central core close to 
the top of the dune. This finding contributes to the understanding of sedimentary structures found 
in nature and might be helpful to improve existing dating methods for desert dunes. 
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If a fluid streams over an extended area of sand, the 
grains will self-organize by forming complex granular 
structures like ripples or dunes [l| . The dynamics of these 
systems is determined by the individual fate of the par- 
ticles 0| . In general, agitated granular matter is known 
to show de-mixing whenever particles differ in size or 
density and indeed size segregation is a well known 
feature of ripples and dunes [4j-ll6j| , as reviewed by Klein- 
hans (l7| . 

Already in 1993 Anderson & Bunas [8j demonstated 
the effect of size segregation in a migrating dune by mod- 
elling the trajectories of large and small particles with a 
stochastic cellular automaton. This kind of particle sort- 
ing observed in natural dunes is known to geomorpholo- 
gists as 'reverse grading', where larger particles are more 
likely to be found in the upper parts of a dune. These 
coarsened crests occur only for fully developed dunes, 
when their lee slopes are shadowed from impacting par- 
ticles. Then, the trajectory of an individual grain behind 
the dune's crest depends strongly on the wind velocity 
field and its mass. Larger, and therefore heavier, grains 
travel in small jumps of the order of a few grain diame- 
ters, while smaller particles are able to leap over the crest 
far down into the shadow zone. Consequently, smaller 
particles end up being buried by larger ones at the toe of 
a dune. 

In addition, Makse [To| showed that size segregation 
due to different hopping lengths in the wake of a dune 
competes with a so-called shape segregation during trans- 
port and rolling of particles with different roughness 
along the dune's surface. These two different segrega- 
tion mechanisms in aeolian sand ripples lead to a so- 
called 'inverse grading climbing ripple lamination' and 
'cross-stratification patterns', respectively. The origin of 
the cross-stratification structures reminds to recent ex- 
periments of avalanche segregation of mixtures of large 
faceted grains and small rounded grains poured in a ver- 
tical Hele-Shaw cell 18, HJ. 



or shape but in density the question arises whether a sim- 
ilar de-mixing of heavy and light grains can be found like 
in other granular systems [3]. Up to now this process of 
density segregation during dune migration has never been 
observed experimentally. By using a bi-dense mixture 
of equally sized particles we demonstrate the counter- 
intuitive effect that the heavier particles accumulate in 
the central core close to the top of the dune. This re- 
sult provides a clue towards an advanced understanding 
of the sedimentary structures found in nature 2ol-22|. 
It might be helpful to improve existing dating methods 
for desert dunes [H, Q and can illuminate the origin of 
placers of minerals in various sediments [25 ]. 

Our experimental setup mainly consists of a narrow 
water flow channel, which allows the experimental inves- 
tigation of downsized slices of transverse dunes [26|, [27J . 
The width of the channel amounts to 6 mm and the 
height is 6 cm. We measure the fixed unidirectional flow 
velocity of 0.60 m/s with an ultrasonic Doppler velocime- 
ter. The corresponding Reynolds number of Re = 36000 
is calculated with the height of the channel. After the 
flume is filled with distilled water, two different types of 
beads are poured into the channel. One type is made 
of glass with a density of p g = 2.5 g/cm 3 . The other 




If, however, the particles inside a dune differ not in size 



FIG. 1: Temporal evolution of a bi-dense dune: five time- 
sequenced snapshots showing side views of a developing dune 
slice composed of glass beads, appearing grey, and white ce- 
ramic beads. The direction of the steady water flow is from 
left to right. 
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FIG. 2: (color). Initial preparation and final segregated state 
for three different starting configurations: (a) heavy parti- 
cles on bottom, (b) heavy particles on top, and (c) mixed. 
Panel (d) shows the mass distribution within the common 
steady-state dune slice. The false color code corresponds to 
the density of the particles, i.e. red represents heavy and blue 
light material. A median filter smooths the rendering of the 
beads in (a) to (c). The height of each panel corresponds to 
4.3 mm in the experiment. 



ones are white ceramic beads with a density of p c = 6.0 
g/cm 3 . Both particle types have the same mean diame- 
ter of 0.5 mm in order to exclude size segregation in the 
form of inverse grading Q • A camera records side views 
of the dunes as shown in Fig. [TJ After suitable image 
processing one can clearly distinguish between the black 
background, the glass beads appearing gray, and white 
ceramic beads. 

For the experiment shown in Fig. [T]we choose a fifty- 
fifty by volume mixture of glass and ceramic beads. The 
ceramic beads are poured into the channel before the 
glass beads are filled in. The resulting starting configu- 
ration is shown in Fig. [I] at t = s. After one minute of 
continuous and unidirectional water flow has elapsed this 
initial configuration is inverted. The lighter glass beads 
have settled at the bottom of the dune while the heavier 
ceramic ones appear at the top. Moreover, the experi- 
ment reveals the nature of this paradox effect, which is 
due to three mechanisms, i) The first one is similar to an 
effect known from agriculture, namely, the separation of 
the wheat from the chaff. In our case the lighter parti- 
cles are dragged by the water flow on longer flights than 
the heavier ones, ii) Further downstream the glass beads 
are caught in the wake region due to the characteristic 
recirculating flow [28] and form an elongated nose as in- 
dicated in Fig. [TJ iii) They provide the carpet for the 
dune to roll over. This explains why the lighter particles 
always end up at the bottom. 

To test the universality of the observed segregation ef- 
fect we prepared two other initial starting configurations, 




FIG. 3: (color). Final segregated states for three bi-dense 
mixtures with different volume fractions ip of ceramic beads. 
The color code corresponds to the one in Fig. [2] The height 
of each panel corresponds to 4.3 mm in the experiment. 



namely, the inverse configuration with the lighter parti- 
cles at the bottom (see Fig.[2fb)), and a mixed configura- 
tion. The mixed configuration is the attempt to achieve 
a homogeneous mixture of the two types of beads, which 
succeeds in a beaker, but fails during filling this melange 
into the channel. The different densities lead to a configu- 
ration as shown in Fig.[2fc), where the heavier beads sed- 
iment preferably in the center. The overall mass at each 
of the three experiments remains constant and amounts 
to m = 22.1 g, as the volume fraction of ceramic beads 
is <p = 50 %. 

Making use of the fact that the three initial configura- 
tions achieve one common steady state we averaged the 
images during the last five seconds of each measurement. 
The result is shown in Fig.[2jd) and represents the attrac- 
tor. The color code indicates the density distribution of 
the particles, which is similar to the one shown in Fig. [T] 
at t = 60 s. 

Another option to test the robustness of the observed 
segregation effect is the change of the volume fraction ip 
of dense ceramic beads. This means that the total size of 
the dune is kept constant, while the percentage of heavy 
beads is varied. Figure [3] shows that the heavy beads 
always end up on top of the dune. However, the contours 
of the dunes differ slightly. At ip = 25 % the downstream 
nose is more pronounced and the slipface formed by the 
ceramic beads is smaller. 

To characterize the dynamical behavior Fig. |4] shows 
the dune velocity v as a function of its mass. It is known 
that in the scaling limit the velocity will scale with the 
square root of the mass, if the dunes consist of only one 
particle species [29|. In particular, it was demonstrated 
that this scaling holds for at least one and a half decades 
and for three different Reynolds numbers 26, 27}. The 
data shown in Fig. |4] for pure glass and pure ceramic 
dunes, respectively, are in agreement with that fact as 
indicated by the dashed black lines, which are fits to 
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FIG. 4: (color). Linear (a) and double logarithmic plot (b) of 
the migration velocity versus dune mass for dunes composed 
entirely of glass beads (squares) or ceramic beads (circles). 
For bi-dense dunes (full discs) the volume is kept constant, but 
the volume fraction ip of ceramic beads is varied, which leads 
to different dune masses. Dashed black lines are power law 
fits with an exponent of -1/2, whereas solid red lines indicate 
power law fits with an exponent of -3/2. In the inset the dune 
velocity v of bi-dense dunes is replotted as a function of the 
mean particle density /Wan- 

However, for dunes consisting of more particle species 
of different density we expect a scaling according to 
v ~ Pmean? with the mean density of the mixture given 
as Pmean = Pg + ■ (pc — Pg), see inset of Fig. HJ This 
exponent can be explained with the finding of Bagnold 
[l| that the migration velocity of a steady-state dune is 
inversely proportional to the density p me an of its sand. 
Combining both ideas we get v ~ Pmean ' ^ _0 ' 5 ? which 
finally leads to 

v-m- 1 - 5 or v~p-l* n (1) 

respectively, under the restriction of a constant volume. 
The corresponding fits are shown as red lines in Fig. HJ 

In addition, Fig. 5 shows the results of an analysis of 
the same data as displayed in Fig. 4. Here, the dune ve- 
locities for a fixed mass m = 22.1 g (see Fig. 1) have been 
derived from the prefactors of the scaling laws for the mi- 
gration velocity v vs. dune mass m. It can be clearly seen 
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FIG. 5: Migration velocity v versus normalized inverse av- 
erage density Pgiass//Wan for dunes composed of glass beads 
and ceramic beads. Here the mass m = 22.1 g is kept con- 
stant, but the volume fraction ip of ceramic beads is varied, 
which leads to different dune volumes. The solid line is a 
linear fit to the data. 

that, for a constant mass, the migration velocity scales, 
as suggested by Bagnold [l| in Chapter 13, inversely pro- 
portional to the average density of the dune's particles 
and thus linearly with the volume of the dune. 

In conclusion, our new measurement technique allows 
to study segregation phenomena both outside and inside 
migrating dunes with high spatial and temporal resolu- 
tion. This method can even be used for tracking all in- 
dividual particles as demonstrated in Ref. Q, which re- 
vealed the inner dynamics of a barchan dune slice. Our 
experiment presents the counter-intuitive effect that in 
bi-dense packings the heavier grains accumulate at the 
top of migrating dunes while lighter particles are buried 
at the bottom. If the particles' densities were distributed 
continuously results comparable to the inverse grading of 
polydisperse packings would be expected. 

As a side effect we show that the migration velocity of 
bi-dense dunes scales with the mean density of the grains 
as a power law function with an exponent of -3/2. This 
experimental observation raises the question whether the 
migration velocity only depends on p m ean or whether the 
distribution of densities also plays a role. 

Another issue is the case of transverse dunes. For tack- 
ling this task new experiments with completely covered 
floor are needed. There, the dynamics of a rippled surface 
developing from an initially flattened bottom has been 
studied to a great extent using a single particle species. 
For bi-dense packings it is expected that similar segrega- 
tion phenomena are observed as for single barchan dunes, 
with the additional feature, that lighter particles leaving 
the lee side of dune will eventually also cover the wind- 
ward side of its downstream neighbor. 

For further examining the underlying mechanisms of 



the segregation phenomena observed so far one has to 
distinguish between size, shape and density segregation 
by examining differences in suspended flight behavior be- 
tween larger /smaller particles and denser/less dense par- 
ticles. Based on our experimental findings, it might be 
speculated that the particle's mass is the most relevant 
parameter. Here the properties of the driving fluid also 
play a role: e.g. for lower Reynolds numbers different 
mechanisms of entrainment of particles into the fluid and 
their successive redeposition on the dune's surface might 
result in a different outcome of the segregation patterns. 

Loosely speaking, these new insights into the sedimen- 
tology of dunes composed of different types of sand has 
the implication, that in a ripple or dune mixed of gold 
and sand, the gold nuggets are likely to be found at the 
top, close to the surface at the crest. 
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